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I discuss in this paper the continuum limit of integrable spin chains based on the 
superalgebras sl{N/K). The general conclusion is that, with the full "supersymmetry" , 
none of these models is relativistic. When the supersymmetry is broken by the generator of 
the sub u(l), Gross Neveu models of various types are obtained. For instance, in the case 
of sl(N/K) with a typical fermionic representation on every site, the continuum limit is the 
GN model with N colors and K flavors. In the case of sl(N/l) and atypical representations 
of spin j, a close cousin of the GN model with N colors and j flavors with flavor anisotropy 
is obtained. The Dynkin parameter associated with the fermionic root, while providing 
solutions of the Yang Baxter equation with a continuous parameter, thus does not give 
rise to any new physics in the field theory limit. 

These features generalize to the case where an impurity is embedded in the system. 
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1. Introduction 



The continuum (field theory) limit of quantum spin chains with ordinary symmetries, 
whether integrable or not, is generally well understood, and described by field theories 
whose symmetries closely match the ones of the underlying lattice model. The most 
striking example of this phenomenon is furnished by the case of sl(N) integrable spin 
chains, whose continuum limits are SU(N) Wess Zumino models, with a level that depends 
on the representation used to build up the chain 0. Quantum group deformations of the 
lattice symmetries are also known to give rise smoothly to similar deformations in the field 
theory and the associated scattering matrices. 

In contrast, the continuum limit of spin chains based on super algebras is rather poorly 
understood. This is unfortunate, since the question is related to physical problems of the 
highest interest, in particular in the context of disordered systems [|J, @, ||, and maybe 
of N — 2 supersymmetry [Kj. One thing that seems clear, is that, in the integrable case, 
none of these continuum limits have to do with the corresponding Wess Zumino mod- 
els on supergroups: this is only expected, since these WZW models can present very 
pathological non unitary properties 0,0: if not, what, then, are these continuum limits? 
The same question arises after quantum group deformation. Here, some preliminary re- 
sults have indicated a very rich structure: it was indeed shown in H that the continuum 
limit of a particular osp(2/2) q model || coincided with the continuum limit of the well 
known Bukhvostov Lipatov model [|10|], providing the first example of an integrable "dou- 
ble sine-Gordon" model O], |E|. The case of sl(N/K) q models based on fundamental 



representations has also been studied in some details [|13|], enough to show that they have 
a relativistic limit, but far from providing a complete identification of the latter. It is thus 
pretty clear that a bunch of solvable field theories of the highest interest are lurking behind 
super spin chains, and this paper is a first step at clarifying the situation. 

Putting the question of field theory aside, integrable lattice models based on super 
algebras have a rather long history, starting with the t-J model, which corresponds to 
si (2/1) and the fundamental representation. Generalized t-J models, based on sl(N/l) 
and still the fundamental representation have also been studied in the context of strongly 
interacting electrons, and - although maybe the algebraic origin of the models was not so 
clear - in the study of quantum impurity problems with the degenerate Anderson model. 
Following developments in superconductivity, models based on more complex algebras or 
representations have been considered: for instance, the model based on the algebra sZ (2/1) 
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and typical four dimensional representations was introduced in [TJ], while the model based 
on si (2/2) and the fundamental (which is also typical) was introduced in [15 . 



The paper is organized as follows. In sections 2, 3 and 4, 1 introduce the general Bethe 
ansatz for chains based on typical fermionic representations. The ground state and physical 
excitations are studied in section 6, while in section 7, I determine the S matrix and the 
mapping onto the Gross-Neveu model. Sections 7 and 8 are devoted to special cases, in 
particular atypical fermionic representations and "bosonic" representations. The whole 
study is extended in section 9 to the case of chains with impurities. A few conclusions are 
gathered in section 10. 

Before starting, I would like to stress that, although the continuum limit of these 
models has not been systematically studied before, the following has a some overlap with 
known partial results from [16| and [17 . 



2. The Bethe equations 

Recall that sl(N/K) has N -\- K — 1 Cartan generators, the first N — 1 belonging to 
sl(N), the last K — 1 to sl(K), the special generator Hn being associated with the odd 
root. The Dynkin diagram decomposes into sl(N) and sl(K) parts, connected by the odd 
root: 

ai a2 &n cln+k-i 

o- -o- - - <g>- - - o- -o 



I wish to consider first integrable hamiltonians with a fermionic representation of sl(N/K) 
on every site. This is a priori the most interesting case, since these representations exhibit, 
in the typical case, a continuous parameter - a feature that is absent in models based on 
ordinary algebras. One might hope that this parameter describes some interesting new 
physics - maybe giving rise to a multiparameter integrable quantum field theory. As we 
will see shortly, the detailed study of the exact solution does not support that expectation 
unfortunately. 

The Dynkin parameters of what I call (a bit incorrectly) fermionic representations are 
(0, . . . , 0, t, 0, . . . , 0), t being a real number, which I assume positive in what follows: the 
case t < would follow simply by exchanging N and K. For t generic, this representation 
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is typical, with dimension 2 NK , and vanishing super dimension. Atypical cases correspond 
to t integer, -{K - 1) < t < N - 1. 

I call the roots of the Bethe equations fiN-i, ■ ■ ■ > A^ij A^o? Ai, . . . , A^-i, and introduce 
the function 

, s v + it/2 , . 

= imp- t 2 - 1 ) 

The Bethe equations read then 

1 =JJe 2 (/UAr_i - fx' N _ 1 )e-i((x N -i - /ujv-2) 

1 = JJe_i(>p - ^ p _i)e 2 (/i p - y! v )e-x{[i v - Hp+i), p = N - 2, . . . , 1 
e^Oo) = JJei(>o - m)e-i(fio - Ai) 

1 = JJe_i(Ai - A 2 )e 2 (Ai - A' 1 )e_i(Ai - /x ) 

1 =JJe_i(A p - Ap_i)e 2 (A p - A p )e_i(A p - A p+ i), p = 2, . . .,K - 2 
1 =J^e_x(Air_i - Ak_ 2 )c 2 (Ak_i - A'^.J, 



(2.2) 



where as usual ||18|| , the pattern of e labels reproduces the Cartan matrix: the two salient 
features are the absence of /aq, (jlq coupling, and the opposite couplings of fio to \i\ and Ai 
respectively. The notation is obvious but implicit. For instance in the first equation the 
product is taken over all Bethe roots and all Bethe roots li' n _ 1 different from Hn-i- 

To get some intuition about these equations, one can think of the dimension of the 
typical representation 2 NK as the number of possible ways of putting fermions with N 
colors and K flavors on a given site of the chain. For instance in the case of sZ(2/l), the 
four states can be interpreted as empty, one fermion with spin up or down, and finally a 
pair of fermions. The parameter iV Mo can be interpreted as the number of fermions; the 
numbers of fermions with a given color are then given by iV Mo — iV Ml , iV Ml — iV M2 , . . . , N^ N _ 2 — 
AT AlJV _ 1 , AT AlJV _ 1 , and the numbers of fermions with a given flavor by iV Mo — N\ 1 ,N\ 1 — 
N\ 2 , . . . , N\ K _ 2 — N\ K _ 1 , N\ K1 . Dynkin parameters can easily be deduced from this and 
the knowledge of the Cartan generators in the typical representation. 

The energy takes the form 

where I have put a chemical potential for the number of fermionic Bethe roots, and e = ±1. 
Explicit expressions for the hamiltonians themselves can be found in the references below 
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for some special cases; they are, of course, quite intricate, except for the simplest values 
of N and K. 

I am not aware of a general derivation of equations ( [2.2| ), even though it is presumably 
possible using the general techniques developed in |T9[], and the form is very natural from 
algebraic considerations |18[]. A number of particular cases have already been studied 
however; besides s/(2/l) |2(J and sZ(2/2), recall that the fermionic representations we 
are considering can become atypical for special values of the parameter t. In the case of 
sl(N/l), the value t = 1 corresponds in fact to the fundamental representation, and the 
model we are interested in coincides then with the su(N) t-J model which was extensively 
studied by Schlottmann [^TJ. There also exists by now a huge literature of quantum 
deformation of super groups and various considerations about graded inverse scattering 
method, with motivations ranging from properties of electronic materials to knot theory. 

The solutions of the Bethe equations are as follows. Consider first the fermionic Bethe 
roots /j,q. Because there is no (/io, Ho) coupling on the right hand side of the Bethe equations 
(the corresponding element of the Cartan matrix vanishes), the usual string solutions, well 
known for ordinary algebras, are not possible. However because the coupling between no 
and Hi has a sign opposite to the one for ordinary algebras, it is possible to compensate for 
the growth or decay of the left hand side of the Bethe equations when has an imaginary 
part by having complexes of "strings over strings". Such complexes were probably first 
introduced by Takahashi |22[ in his study of one dimensional fermions interacting with 
an attractive delta function potential; they have been widely used since, in particular by 
Schlottmann in his study of models based on the fundamental of sl(N/l). They are of the 
type 

^ = if" 1 + (-(p - l)i/2, . . . , (p - l)i/2) 

ft q = yP~ x + (-(P - 1 - q)i/2, . . . , (p - 1 - q)i/2), q = 1, . . . ,p - 2 (2.4) 
Hp-i = (J? -1 , 

for p = 1, . . . , N. Hence, for (jlq there are string solutions of length smaller or equal to N. 
Of course the patterns ( |2.4| ) are obeyed only in the large L limit. For finite L the solutions 
of the Bethe equations differ from ( |2.4[ ) by exponentially small amounts. One has to use 
and eliminate these small deviations to rewrite Bethe equations involving the complexes. 
As for the /j, p , p > roots that are not involved in such complexes and the A p roots, they 
are determined by the same arguments as for sl(N) and sl(K) respectively ie they can be 
strings of any possible length. Observe that the role of sl(N) and sl(K) are exchanged if 
t is negative; once again, in the following, I shall assume that t > 0. 
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By taking the logarithm of the Bethe equations and differentiating we get a system 
of integral equations. Let us introduce notations for densities. I call p p (p = l,iV) the 
density per unit length of real centers of fxo strings of length p (|2.4| )(ie the density of 
in ( |2.4p ). The density of real centers of / strings of \x v roots that are not in one of the 
complexes ( j2.4| ) we call Op . The density of real centers of / strings of X p solutions we 
call Tp \ I will usually reserve the labels p, q for the colors of roots and Z, m for the types 
of strings solutions. We also use the labels p, q for the complexes of strings over strings 
because they behave in many ways like new roots colors. 

I define the Fourier transform as 

/(*) = J dve^f(v), /(i/) = i- J dxe-^f(x), (2.5) 

and introduce the following notation 

a t (u) = — 4- In \e t (u)] = - t — Tl (2.6) 

K ' 2ndv 1 K n 2n u 2 + f 

with 

d t (x)=e- t ^/ 2 . (2.7) 
I also define for r, s integers (all these notations are rather standard) 

G rs = a r _|_ s _2 + a r _|_ s _4 + . . . a r _ s , r > s; G rs = G sr , r < s, (2-8) 

(with ao = 0) and 

A rs = (G rs + 5 rs )*(l + a 2 ). (2.9) 
where * denotes convolution. Their Fourier transforms are 

6r, = ^^e-^ft-6„, r>s, (2.10) 
smh.{x/2) 

and 

j 2cOSh(x/2) -HtI/2 ^ rn n\ 

A " = • u/ L sinh xs/2 e r|a;|/2 , r > s, 2.11 
smh(x/2) 

I also introduce the kernel 

/ x tv/2 , . 1 , 

= r , = : --, 2.12 

v ' coshyrz/ v 7 2coshx/2' v ; 

and define 

Clj-s ~~ S ~k ^4. rS; (2.13) 
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with 

A , . sinh(xs/2) _ r \ x \/ 2 ^ , n . , N 

M*) = ■ A L e r|x|/2 , r > s. 2.14 
smh(x/2) 

I can now write the continuum version of the Bethe equations. Introducing the symbol 

Gt+i, p which is defined by a formula similar to ( |2.8| ) even when t is not integer (which is 
usually the case) 

Gt+i, P = at+p-i + at+p-3 + ■ ■ ■ CLt-p+i, (2-15) 

(so, for instance, Gi, v = a p ) we have, for the fermionic root, a set of iV equations, one for 
each string 

N 

G t +i, P = p P + pp + ^Gpq* pq + ^ai* a p l) -^2,a p i*T{\ p=l,...,N. (2.16) 

q=l 1>1 1>1 

For p = N recall that there is no density a^j so the corresponding term has to be suppressed 
from the equation. For the N — 1 roots of sl(N) we have an infinity of equations, one for 
each type of string 

N-l 

^*P P = ^ ) +E A '™*E^* (T i ra) ' p=l,...,W-l, 1>1 (2.17) 

m>l <Z=1 

where 

C pq (u) = 5{u)5 pq - S(v)8p 7 q-1 - s(u)S p>q+ l. ( 2 -18) 

Finally for the K — 1 roots of sl(K) roots we have 

N K-l 

S P i^a lq *p q = f® + J2 A im^J2 C Pi* T t^ P=h---,K-1. (2.19) 



3. The thermodynamic equations 

We now write the thermodynamic equations. The energy term is 

N 

E = J2 (eG t+hp +pA)p p , (3.1) 
p=i J 

where for notational simplicity, we have not written the variables that are integrated over 
(the real centers of the fermionic strings). The chemical potential breaks the "supersym- 
metry", leaving as a symmetry the bosonic part sl(N) ® sl(K) <S> u(l). 
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We introduce pseudo energies e p , n p l \ (jjp defined by 

Pp/Pp = exp(-e p /T), p = 1, . . . , N 
a (D/aW = exp(-4')/T), p = 1, . . . , N - 1, / = 1, . . . , oo (3.2) 
r«/f« = expC-C^/T), p = 1, . . . , K - 1, / = 1, . . . , oo. 

The minimization of the free energy leads to the system of thermodynamic Bethe ansatz 
(TBA) equations: 

N 

. , 4/V'/- 

q=l 1>1 

-J> p *ln (l + e-^/ T ) 
l>i 

and 

N-l 

= -Tin (l + e4°/Tj + £ Ami * J- C qp *Tln (l + e~<' '/ T ) +a t *T]n (l + e"^ T ) , 

(3.4) 



=pA - G t +i,p -e p + ^G pq *T\n(l + e- £ ' /T ) - ^ aj ★ Tin (l + e"^ 



(3.3) 



m>l q=l 

and 



= - Tin (l + etfVr) + J] ^ * E ^ * Tln ( X + e- C ^ /T ) 

m>l 

JV 

-V^a^rin ^i + e -^/ T ). 



^ ^ (3-5) 



The free energy reads then 

(3.6) 



f =4/^ G «' k ( 1 + e ""'1 



4. Large temperature entropy 

Before going any further, it is useful to check the completeness of the solutions by 
studying the large temperature entropy. The general case is a bit heavy, so I will sim- 
ply discuss some particular examples here. Let us start with si (2/1), and introduce the 
quantities 

xi =e K i I , / = l,,...,oo 

Vi =e £l/T (4-1) 



V2 =e^ T . 
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In the large temperature limit, these go to constants, which are solution of the system 

Sn/2 

) 



Xl = [(l + x,_i)(l + x l+1 )] 1/2 + - 

V Vi 



2/i 



1)2 



1 + 


1 






2/2 




1 


1 + 






2/2 



(4.2) 



The solution of this system is 



i + 3 - 

2 



1 

2/1 



1, 2/1 = 2/2 = 3. 
5 



Meanwhile, at large temperature, one has, for t ^ 1 



-T 



2 In 1 + 



1 

2/2 



+ In 1 + 



2/1 



-Tin 4, 



(4.3) 



(4.4) 



in agreement with the dimension of the typical representations of sZ(2/l), d = 4. If, 
however, t = 1, one finds 



T « -T 



In 1 + 



1 

2/2 



+ In 1 + 



1 

2/1 



-Tin 3, 



(4.5) 



in agreement with the dimension of the fermionic atypical representations of s/(2/l), d = 3. 

wehave performed the same exercise for sZ(3/l). It is a bit more laborious, so wewill 
only give the final result here. Introducing 2/i = 2/i = e £i / T , wefound y\ = |, z/2 " 
r/3 = 7. The large temperature free energy for typical representations is then 



17' 



1 



2/3 



1 



3 In 1 + — +2 In 1 + — + In 1+ — 



2/2 



1 



2/1 



-T In 8. 



There are now two types of atypical representations. If t = 2, one has 



T ! 

while if t = 1, 



-T 



2 In ( 1 + — ] + 2 In ( 1 + — 

2/3 / V 2/2 



F w -T 



1 \ 



1 \ 



-Mi) ( 1 + — J 

2/1/ 



1 \ 



-Tin 7, 



In [ 1 + — + In 1+ — + In 1 + — 

2/3/ V. 2/2/ V. 2/1/ 



-Tin 4. 



(4.6) 



(4.7) 



All of these coincide with known results of s/(3/l) representation theory. The general 
relation between the sl(N/K) TBA and representation theory seems quite interesting, but 



I won't comment any more on it here [ 23] . 
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5. The ground state and physical excitations 

In this paragraph, I will restrict to the "generic case" t > N. Some special cases are 
studied further below. I will also concentrate on the case e = — 1, and comment briefly on 
the case e = 1 - which happens to be quite similar - at the end 

5.1. Equations as T — > 

As T — > we find the system 

N 

= pA- G t +i, p - e p -^2G pq *e~ +^2ai* K p l) ~ +^ai p -k Q l '~ , (5. 1) 

q=l l>\ 1>1 

N-l 

o = - A mi* E c ip * ~ a l*tp' (5-2) 

m>l <j=l 

K-l N 

= -Cp )+ - E A ^ * E c ™ * ^ m) ~ + 5 pi E °«« * v ( 5 - 3 ) 

m>l 9=1 9=1 

where I have introduced as usual the positive and negative parts of the pseudoenergies. 
Whatever the value of A it is easy to see that one has 



and 



and 



together with 



and 



ei = ... = e N _ 1 =0, (5.4) 



K (0± =0, p = 1, . . . , JV - 1, I = 1, . . . , oo. (5.5) 



C«+=0, Z = l,...,oo,; C^"=0, MiV. (5.6) 
The equation (|5.2| ) is then satisfied identically. The equation ( |5.3|) now reads 

K-l 

= -A lN ★ E * C^ - + <Wjvz * (5.7) 

9=1 

Using a^; = s ★ A/vz ( |2.13|) , eq. (|5.7|) can be rewritten 

K-l 

= -EV Cf " + <W * (5.8) 

9=1 
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that is, using the form of C 3 



(AT)- a ^AN) 



cr~-**c. 



AN)- 
S2 



S * 



2 — s -k e N 

4 N) -+ci N) - 



(5.9) 



Sfc-l — * * SK-2 i 



whose solution is 



sinh(X — p)x/2 
sinhi^x/2 



-N- 



(5.10) 



or 



SK-p,K * Cjv, 



(5.11) 



where I defined 

sinhrx/2 
rs sinhsx/2 

Similar results hold for the densities, that is 



(5.12) 



Pi 



Pn-i = 0, 



(5.13) 



together with 



r«=a«=0, p = l,...,7V-l, Z = l,...,oo, 



(5.14) 



and 



fW=0, p = l,...,K-l, / = !,.. .,oo, 



(5.15) 



and 



-«=0, p=l,...,X-l, Z^iV. 



(5.16) 



From (|5.11j ) we get also 



Tp N) = S K -p,K* PN- 



(5.17) 



We thus see that, whatever the value of A, the only non- vanishing particle densities are 
Pn and Tp N \ To proceed further we have to distinguish the cases A = and A > 0. 
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5.2. The case A = 

I consider first the case A = 0. In the ground state one has then, in addition to ( |5.4| ) 



4 = •••4-1=0, (5.18) 

so all the hole densities vanish, leaving the system 

= -G t +i >p - (S pN + G pN ) * + a Np * ([ N) ~ , (5.19) 

Prom ( |5.11[ ) we replace Q by its expression in terms of to get, in terms of Fourier 
transforms, 



Q = sinhpx/2 c _ f|x|/2 _ smhpx/2 c _ (N _ 1)]x]/2 



sinh x/2 sinh x/2 N 

sinhpx/2 _ N \ x \/ 2 sinh(iiT — l)x/2 

sinh x/2 sinhXx/2 

As expected, p disappears and we get 



(5.20) 



? _ _ sinhKx/2 j N _ K _ t - ) \ x y 2 un)- _ sinh(K - p)x/2 jN-K-t)\x\/2 /c 91 \ 
" N ~ sinh x/2 6 ' ip ' ~ sinh x/2 6 ' ^ J 

and therefore 

e N = —a t -N+K,K, C,p N ' > = —at-N+K,K- P , (5.22) 

all other pseudoenergies being zero. In the ground state we therefore have as well, by 
comparing the Bethe equations and the limit T — > of the thermodynamic equations, 

= ^Kx/2 eiN _ K _ t){w{/ ^ = 
sinh x/2 

• UJT \ 10 ' (5 - 23) 

-(AT) Sinn^A — p)X/Z JN-K-t)\x\/2 _(AT) 

^ = sinhx/2 ^ > ^ = 

all other densities being zero. The ground state is thus filled with complexes of iV strings 
over strings (p.4j ) and N strings for all the sl(K) roots X p . 

Excitations are made of holes in the distributions (|5.23|) , with excitation energies 



exactly equal to — and —(p N ^ ■ By a standard argument the momentum is given by 
p = 2-k J (density). Taking inverse fourier transform we see that these excitation energies 
are expressed as sums of terms of the form atiy) while momenta are sum of terms of the 
type ilnet(y). At large values of the bare rapidity v where the gap vanishes we therefore 
have e oc and p oc - ie e oc p 2 . The excitations are therefore non relativistic. Such a 
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dispersion relation is characteristic of quantum ferromagnets. However I have not chosen 
the "wrong sign" of the hamiltonian: identical features are observed for e = 1 (see below). 
There does not seem to be a very physical reasons why the excitations are not relativistic. 
Technically, what happens is that the fermionic Bethe roots having no self coupling, the 
dispersion relation of the associated dressed excitations is almost the same as the one 
of the bare excitations. This can be seen especially clearly in the case of si (1/1) (more 
generally, N = K), where the Bethe equation reduces to et(^o) L = 1, and the energy 
of excitations is = —at- The hamiltonian is the one of a XX chain with a magnetic 
field, H = (&j^cr~ + cr j~ cr j~) ~ ^J2j cr j- After fermionization and Fourier transform, 
it becomes H = 2 cos ka\ak — 2F, F the number of fermions. In that language, the 
ground state is obtained by filling up all modes — n < k < n, and the gapless excitations 
occur near k = 0, where the energy goes like e oc k 2 . 

I thus conclude that, if the supersymmetry sl(N/K) is not broken, integrable lattice 
models based on fermionic representations do not have a relativistic limit (we will see later 
that this is true for other representations as well). This result is of course disappointing, 
and in sharp contrast with the situation for ordinary Lie algebras. To get some non trivial 
results, we do in fact need to break the supersymmetry, as I now demonstrate. 

5.3. The case A > 

Suppose now A > 0. The first difference with the case A = is that e has also a non 
vanishing positive part obeying 

G t +i,N -NA = -e+ - (S + G) * e^, (5.24) 

with the kernel 

~_ smh(N-K)x/2 N]x]/2 

G " sinhKx/2 6 • {b ' 2b) 

In particular when iV = K we have simply (much like in the N = K = 1 case) 

G t+ljN -NA = -e+ -e„, (5.26) 



The function e^v is now negative on a finite interval [-Q, Q}. We have 

-Q 
-Q 
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rQ 

eJfiX) + G(X - v)e^(fi) = -G t+1 , N + NA, A e [-Q, Q], 

J-Q 



(5.27) 



and 

4+ i 

|A»|>Q 



e++ H(X - p)e%(p)dp = —a t _ N+KjK + KA, \X\ > Q (5.28) 

J\u\>Q 



where 1 + H = — ^ . 

l+G 

These equations can be solved perturbatively in the limit Q » 1 using standard 
Wiener-Hopf techniques. At dominant order in this limit, the cut-off Q is related to the 
magnetic field by 



A « —a t - N+K , K (Q) « , Q » 1, (5.29) 



1 est 

— a t _jv+K,K(Q) « ^ 

(this result is exact in the case iV = K). The system is still gapless but the gap now 
vanishes at finite rapidity. For p close to ±Q one has 



e N (p) ^\a' t _ N+K>K (Q)\(\p\-Q), \p\ » Q. (5.30) 

In the presence of a magnetic field there is no simple relation between the ground 
state pseudoenergies and densities. The latter obey, instead of ( |5.27| ), 

rQ 

Piv(A) + / G(X- p)p N (p)dp = G t+1>N , Ag[-Q,Q], (5.31) 
J-o 



and 

Pn 



+ / H(X- p)p N (p)dp = a t -N+K,K, |A|>Q. (5.32) 

J\v\>Q 

In particular these densities are discontinuous at the cutoff Q. For large Q we have ap- 
proximately Pn(Q) ~ Pn(Q) ~ dt-N+K,K(Q)- On the other hand the relation between 
momenta and densities still holds so we get for the momentum of excitations (|5.30|) 

Pn(p) ~ 2na t -N+K,K(Q)(\p\ - Q), \p\ ~ Q (5.33) 

The massless excitations in the branch therefore now are relativistic, with the sound 
velocity 

i K_„ W (Q)L. i Q>>L (5 ^ 34) 



2n a t - N+K)K (Q) nQ 
The part induces also non vanishing e+. We have 

= pA - G t +i,p - e+ - s pA r ★ (5 + G) ★ e^, (5.35) 

where we used 

^ _ a sinh(if - l)x/2 _ sinhpx/2 
pN pN sinhiira;/2 sinhXa;/2 
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Observing that 

^ sinhpx/2 
Gt+1 > p = smhNx/2 Gt+1 > N 



and using (|5.24|) we get 



£p=s pN -kej f , (5.36) 

Therefore, we have now new particle like excitations in the system, with energy e+ . 
Their physical nature is easy to understand. At a given rapidity, when Q >> 1 only the 
tails at ±00 of s p n gives significant contributions because e^- vanishes in [— Q, Q]. At large 
argument the behaviour of s p n is determined by the pole of its Fourier transform nearest 
the real axis, that is x = 2itc/N , and we can approximate 

s pN {y) - I sin (^) e-fcrM/i* v ioo . (5 . 37) 
Expanding close to Q, we get therefore 

e+(//) « ^sin(^) |a;_ Arw (Q)|e- 2 ^/ 7V cosh(27r^/iV), \fi\ « Q, p= 1,..., JV-1 

(5.38) 

The momentum of these excitations is given by p = 2n J p p . From ( |5.35|) we deduce as 
well 

P P = s p n * Pn, (5.39) 

and thus 

p = ^ sin (^) a t - N +K,K(Q)e- 27TQ/N sinh(27r^/iV), |/x| « Q. (5.40) 
Hence the excitations are relativistic once again, with a mass 

m P = ^ sin (^) \a>_ N+K , K (Q)\e- 2nQ/N , P = 1, • • • , - 1, (5.41) 



and the same sound velocity as before (which was quite obvious from (|5.36| ) and ( |5.39| ). 



The ( excitations are also modified due to the existence of the cut-off Q. This is 
easily studied using ( |5.11| ) which still holds. The gap still vanishes at rapidities much 
bigger (in absolute value) than Q, but with a different dependence on rapidities. In this 
limit the behaviour of (Cp N ^ 1 is determined by the tail of the kernel ( |5.12|) sk- p ,k- For 
A > (< 0), only the region close to Q(—Q) contributes, so we find 

(C^)" ^ $ sin (^^) l«t-^^(Q)|e- 2 ^e ±2 -/^, (5.42) 
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The momentum of these excitations is given by p = 2tc j 4 N) . From (EU7D we check that 



these excitations are now left and right moving relativistic massless excitations with a 
sound velocity that is still given by (|5.34| ) and a mass parameter 



m„ = — 77 sm 

TV 2 



K (^^) \<-n+kAQ)\*-** Q ' K > V = l,..;K-l (5.43) 



To conclude this section, I would like to notice that the techniques I used are entirely 
similar to the ones developed by Tsvelik in his study of si (2) chains with a magnetic field 
51 



5.4- The case e = 1 

When e = 1 and A = 0, it is easy to see that the solution of the Bethe equations 
is = Gt+i, p for p = 1, . . . , N — 1, while all other pseudo energies vanish: the ground 
state is empty, but there are hole densities for all the fermionic strings. When A = —B, 
B > is turned on, only e^v acquires a negative part, while the relations e+ = s Pj n * e^- 
and (p N ^ = sk- p ,k*£n hold. Things are thus very much similar to the case e = — 1: 
the difference is that is now positive (instead of negative) in a finite interval , and thus 
it is the €p excitations that are massless, while the (p N ^ excitations are massive. In effect, 
the roles of iV and K are thus exchanged. 



6. S matrices in the scaling limit 

We now discuss the way the various excitations interact in the case e = — 1 and A 
a small positive number (similar results would hold for e = 1 and A = —B, B a small 
positive number, up to the exchange of iV and K). In the following we shall be interested 
in the scaling limit of the lattice model. In this limit the three types of excitations we have 
identified (e+ , e^-, C,p N ^ ) become decoupled since they occur respectively for rapidities /i 
such that « Q, \n\ ~ Q and >> Q. We shall generally write equations where this 
coupling has been neglected by the symbol ~. Densities evaluated in the ground state are 
denoted by |o- 

We now get back to the equations that involve densities (and are magnetic field in- 
dependent), and consider first the physics in the vicinity of \fj,\ « Q. In that region, it 
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turns out that physical densities are hole densities r, so our first task is to invert ( |2.19| ) 
and express instead the densities r in terms of hole densities f . We find 

K-l 

B pl * S *p n = T^ -J2 B Pl*J2 Cnm * f t ] i V, q=l,...,K-l, (6.1) 

g=l m>l 



where 



2 cosh t /2 

B - = sinh^sinhW Shlh [(K - " )X/2] SiahkX/2] ' P ^ <8 ' 2) 



and C lm is defined exactly like in ( |2.18|) , but acting on upper indices. Also, when n> N, 
there is no density p n , and the source term disappears from the equation (|6.1|). The next 
step is then to replace t[ in ( |2. 16| ) by the expression ( |6.1| ) for p = 1. We then use the 
last equation ( |2.16|) for p = N to eliminate p^. Replacing in the equations ( |2.16|) for 
p = 1, . . . , N — 1 leads to 



Gt+i, P — (Gt+i,N * H pN * (1 + H NN ) 1 ) + H pN * (1 + H NN ) 1 * pn — 

N-l 

P P + ( H pi ~ H p n * Hn i + Hnn)' 1 ) *p q + ^2ai* a® 

q = l 1>1 

K-l 

+ {H P n * a Nl * (1 + Hnn)- 1 - a pl ) * ^ B lq * ^ ^ * ^ 

Z>1 g=l m>l 



(6.3) 



where we introduced the kernels 

The last term vanishes for I > N; moreover, in the approximation we are considering, it is 
not possible to make holes in the distributions Tq N ^ , which would cost a very large energy. 



The term r q thus disappears from the equation (|6.3| ), leaving a finite set of r densities. 



The constant term on the left hand side of (|6.3|) vanishes identically too. After simplifying 
the expressions slightly, we thus end up with the system 

N-l K-l 

Pp+Pp ~ s pN *p N \ + 2^ Z pq *p q - 22 s K-q,K^ p ' -^ai-kcr^ 1 ' , p = l,..., N-l (6.5) 

q = l q = l 1>1 



where 



_ sinhrx/2sinh(N-s)x/2 Klxl/2 
° rs sinhKa:/2sinhiVx/2 6 ' S *' rs sr " { ™> 
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This has to be supplemented by the equations 

N-l 

a( p l) +J2 Al ™*J2 C Pi* a i m) =ai *Pp> P = h---,N-1, 1>1, (6.7) 

m>l g=l 

and 

K-l N-l 

t& -Y^B^+Y,^*^ ™ 8 k-*>k*P» P = 1 i N ~ 1 i q = l,K-l. (6.8) 

r=l s = l 

We can identify the source terms s pn * p~n\ with p/2tt thanks to (|5.40|) . These 
equations are then easily shown to coincide with the system obtained in the study of an 
su(N) scattering theory at level K. Recall in particular that there are particles of mass 
m p = m sin(p7r/iV) associated with every fully antisymmetric representation and carrying 
a charge which is a weight of these representations. The S matrix is discussed in |25[| : it 
is the tensor product of an sl(N) level K RSOS S matrix and a sl(N) "vertex" (soliton) 
S matrix. This structure is transparent on ( |6.5|) : the fq are the densities involved in the 
diagonalization of the RSOS part of the scattering matrix, while the a p are those from 
the diagonalization of the vertex part. 

For the excitations at rapidities \p\ >> Q, things are a bit simpler, since the densities 
p and a are totally frozen in that limit. The physics is thus fully described by the equations 

K-l 

SpiaiN * Pn\ ~ s f$ ) + Alm * C pi * r i m) " ( 6 - 9 ) 

m>l g=l 

These equations are similar to the ones one would write for a sl(K) lattice model with the 
fully symmetric representation Nui on every site, in the limit appropriate to study the 
massless right moving excitations |26|| . The scattering theory can then be easily extracted. 
This time one has massless excitations of mass parameter m q = msm(q7r/K), and the S 
matrix is the tensor product of an sl(K) level iV RSOS S matrix and a sl(K) soliton S 
matrix. This massless theory is well known to describe the SU(K) level N WZW theory. 
The idea of describing a conformal field theory by a massless scattering theory has a long 
history going back to ]27J . It has been a subject of intense interest recently in the context 
of quantum impurity problems |28| and the quantum KdV equation [p9|. 

Finally the excitations are completely free in this limit, describing a massless U(l) 
degree of freedom. 
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The different pieces of scattering theory found in this section coincide with the known 
results [3C] for the iV colors, K flavors chiral Gross Neveu model. We thus conclude that, in 
the limit of large Q (that is, infinitesimally small symmetry breaking field), the continuum 
limit of the sl(N/K) quantum spin chain with generic fermionic representations obeying 
t > N — 1 coincides with the chiral Gross Neveu model 

£ = ififfrPjf + 9^L jf ^Rjg^R g ^Lkf, (6.10) 

where j the color index runs from 1 to N and / the flavor index from 1 to K. In terms 
of currents, the interaction reads JlJr + Jl^ri wnere J is the £7(1), chirality carrying 
current, J a are the sl(N) currents, with = ip^f (T a y k ipRkf- The parameter g in ( |6.10j ) 
varies as g oc 1/Q, and the true scaling limit is obtained when q — > oo, that is g — » 0. In 
that limit, using the equation determining e±, and well known considerations on dressed 



charges ]3T[], I found that the U(l) degree of freedom has a radius R = -j==- Except when 



471" 



N = K, this is not the radius that is expected for the action ( 6.10 ), and the latter requires 
an additional JlJr coupling to be correct. As is well known, such a coupling does not 
change any of the physical or integrability properties in a significant way. 



7. Some particular cases 

We concentrated so far on the case of fermionic representations with t > N — 1 (which, 
in particular, are typical). When t < N — 1, things can get quite complicated, due to the 
different possible structures of the source term Gt+i tP in the Bethe equations. 

The simplest situation occurs when e = 1. In that case, it is easy to see that when 
A = 0, the ground state is always given by e+ = G t +i, p (all others being zero), irrespective 
of the value of t (and thus excitations are not relativistic). When the field A = —B is 
turned on, is the only fermionic pesudo energy to acquire a negative part, which also 
gives rise to negative parts Q = sk- p ,k * e N- The situation is thus very similar to the 
case t > N — 1: the e+ excitations are massless and described by an SU(N) level k Wess 
Zumino model, while the Q are massive, and described by an SU(K) level iV massive 
theory; the e^v excitations still correspond to a simple U(l) theory. The only difference 
with the case t > N — 1 is that the Fermi velocities of these excitations are not in general 
equal anymore. To see this, let us restrict to the case K = 1 for simplicity. The equations 
for the fermionic roots are then 

= - P B + G t +i, P - e+ - G vN e~ N , p < N - 1 

= - NB + G t+ljN - e+ - (1 + G NN )e^ 
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The pseudo energy e^v acquires a negative part at large rapidities, while one has 

,+ _ G pN ,+ i n G pN G t+ i,N ,„ ^ 

6 P ~ 1 _i_ n e N + Lr t+l,p : n • 

In contrast with the case t > N — 1 where the second term vanishes, the massless relativistic 
region now corresponds to rapidities much larger than Q (the Fermi rapidity for ejv), where 
the first term of ( |7.2| ) is negligible: the behaviour of e p is thus fully determined by the 
second term: it is independent of e^v and Q, and thus is bound to have a different Fermi 
velocity than the massless U(l) excitations. More careful study of this second term shows 
that it reproduces a set of massless excitations with mass parameters proportional to 
sinp7r/iV, for any value of t. The continuum limit of this model is thus the tensor product 
of a level 1 SU(N) WZW model and a U(l) boson, each with its own sound velocity. 
Notice that it is not necessary to take the Q — > oo limit here (since there are no excitations 
at small rapidity to decouple) and as Q varies, the radius of compactification of the boson 
changes. When Q — > (low density limit in the t-J language to be discussed below), it 
goes to the point R = of the c = 1. The whole theory thus reproduces the well known 
system of free fermions with U(l) x SU(N) symmetry (see eg [[32]). When instead Q — > oo 
(corresponding to the limit of half filling in the t-J model), the radius goes to R = \J~^, 
the same value obtained when bosonizing a single free fermion (which corresponds to the 
case N = 1). This value arises also when considering the U — > oo limit of the SU(N) 
Hubbard model |33[]. 

Indeed, the case t = 1, K = 1 is nothing but the so called SU(N) t-J model ]2~T| ]. 
It is interesting to notice here that the Bethe equations we started with coincide with 
those of Schlottmann, and therefore to a point of view where there is an atypical fermionic 
representation on every site. As is well known, other Bethe ansatze can be written for this 
model, as was done first by Sutherland [33]: they correspond to a point of view where there 
is a fundamental representation on every site fl35|. The continuum limit of the t-J model 
has been worked out in [36]: the results are identical to what we found here, although they 
are not formulated in terms of massless scattering, but rather using conformal dimensions 
and dressed charge matrices. 

The case e = — 1 is considerably more involved, and it is not clear where the ground 
state lies as t varies ||19j . An exception is t = N — 1, where, when A = 0, the ground state 
is obtained by filling up the sea of length N — 1 fermionic roots, e^_ 1 = — ^nh^/^ ' e ~ Kx ^ 2 i 
and excitations are, as usual, not relativistic. When a field is added, e^-i acquires a 
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positive part, together with the other e's which obey e+ = sin h(jv- X i^/2 e iv-i ■ The relation 
X ^ = sk- p ,k * e AT-i holds, too. It follows that both (p N ^ and e+ give rise to 
relativistic massless (resp. massive) excitations. In addition however, one has Kj^-i = 
—a; * £jv_i , giving rise to non relativistic excitations: from the field theory point of view, 
this is thus a not very interesting situation. I suspect similar conclusions hold for other 
values of t. 



8. Other representations 

In the K = 1 case, the choice t = 1 corresponds, in fact, to putting a fundamental 
representation on every site. It is interesting to study more generally the case where the 
chain is built up using representations with Dynkin parameters 

j 

o- -o- - - 0- - - o- -o 



The equations then look as in (|2.2| ), except for the source terms: the left hand side 
for the /io root is now equal to one, while the left hand side for the Pn-i root is where 
j is the (integer) Dynkin parameter, j = a\. The su{N) t-J model corresponds to K = 1, 
j = 1 already discussed above. In the presence of a chemical potential, the energy reads 

E = -J2 -r^—z + A E 1 - («) 

Ajjv-i H'N-l 4 MJV-i 

The solutions to the Bethe equations are simpler than in the fermionic case: there are 
the usual strings for every bosonic root, while the no for the fermionic root are all real. 
With the same notations as before (setting pi = p) we now have, going to the continuum 
version 

G j+1> i5 p>N -i = 4° + Yl Alm * C pi * a i m) ~ S P iai * P> ( 8 - 2 ) 

m>l <J=1 

together with 

= P + P-J2 ai * a i l) + J2 ai * T i ] > P=h---,N-1, (8.3) 

i>i i>i 
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and 

K-l 

S P iai*p = fW + Y, a i™*Ys C Pi* T t ) > P = l,---,K-1. (8.4) 

m>l <?=1 

To discuss what is going on, let us consider as an example the case N = K = 2. The 
TBA equations at T = are, setting = k® , C? = C (0 , 

G j+1>l = - -J2 A ™i* « (m) " + *i* e" 

o = - c (l)+ -J2 Ai ™* + a i* e " ■ 

This system can be transformed into 

K (0 = s * + «('+!)+) + (5 ilS * e~ - 

e=A + ^a^K ( ' ) - -^ ai *C (0 ~ (8.6) 
C (0 =^(C ( '- 1)+ + C (m)+ )+fe^e". 

In the ground state, all positive parts of the pseudo energies vanish. It follows that k' 1 ' = 
s * e~ and k^' = —s, C^ 1 ) - = s* e~ . These can be put back in the equation for e, which 
reads in that case (it is particularly simple because N = K here) 

e = A - s*a,j. (8.7) 

When A vanishes, the e excitations are non relativistic. When A is positive, and provided 
it is not too big, e has a positive and a negative part; the corresponding u(l) excitations 
are then massless; calling Q the Fermi rapidity for the e excitations, their sound velocity is 
v = ^ ^ p (Q)^ • From ( |8.6| ), it also follows that = —s, = s*e~ and = s*e~. 

The k excitations are thus also massless and relativistic at rapidities much larger than Q: 
notice however that they have different Fermi velocities: the one for is independent 

of Q (it turns out to be Vj = tc in our conventions) , while the one for k^ 1 ' does depend 
on it, and reads v\ = ^WF^T^ > /U — > oo. As for the ( excitations, they are also massless, 
with a similar sound velocity. 

Much like in the previous section, and in contrast with the case of generic fermionic 
representations (with t > N — 1), relativistic invariance does not require Q to be large, and 
thus is obtained for an entire range of values of A. Here, this means that the associated 
scattering theory has a continuous parameter A, whose meaning we now partly elucidate. 
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To do so, we observe that the equations for the densities are identical to a decoupled system 
for a pair of sl(2) spin chains, the first sl(2) chain having a source term on the first and j th 
node, the second si (2) chain on the first node. The second system has thus the su{2) level 
1 WZW model as continuum limit. As for the first, it is similar to a general class of lattice 
models with mixtures of several spins. These models were discussed in details in [37J and 



38fl . In [5S|, the sound velocities for all excitations were the same, and the continuum 



limit was described as the tensor product of an SU (2) level j WZW model (with c = t-^) 
and of an SU{2) minimal coset model with c = 1 — ^. +1 ^. +2 ^ . In the limit of small Q, 
this result essentially still holds, but now the two types of excitations each have different 



Fermi velocities, Vj and V\ respectively. The scattering theory is as in [38], and the field 
theory can be related with a 2 colors j flavors Gross Neveu model with flavor anisotropy. 
Away from that limit, things are more complicated: the contribution to the free energy 
from the k degrees of freedom can be written as f K ~ — ^f— (jf: + f^Jj arid although the 
sum Cj + c\ stays the same (it is controlled by the overall shape of the TBA diagram), the 
individual values of these two parameters evolve with A. In the limit where Q — > oo, one 

3 ( ' 1 1 

finds that c\ = 1 while Cj = . One should not however think that this model always 

decomposes into the sum of two independent CFTs with different sound velocities, and 
central charges Cj and c\, except for very small and very large Q. To get an idea of what 
happens, I will restrict to the case j = 2, which was also partly treated in [17]. Consider 



therefore a theory made initially of an Ising model and a level 2 WZW model, which is 
the right description of the system at small A. This theory can be written in terms of 4 



species of Majorana fermions Xh i = 0, 1, 2, 3. Suppose now one adds, as suggested in [|T7f] , 
a four fermion coupling in each chiral sector. The hamiltonian say for the right movers 
reads then 

3 

n = -ia^2 XidxXi ~ ibxodxXo + 0:0X1X2X3- 
i=i 

To handle this model, the best is to bosonize, representing the fermions as xo ^ cos 4>ii Xi (X 
sin (pi |3~9fl (and similarly for X2, X3 m terms of another boson 02 ). The hamiltonian reads 
then, schematically, 



n = a 



(d0i) z + 2 (dfcT + iao&fa +B (d0i) 2 -ia d 2 - 



n 



where d = d x . For C = 0, and a choice of «o corresponding to the bosonized Ising stress 
energy tensor, the A term is a c = 3/2 theory, the B term a c = 1/2 theory. Moreover, 
these two theories are independent (the short distance expansion of their two stress energy 
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tensors is regular). For a = 0, and A = £?, the theory decomposes in two independent 
bosons 4>i ± 4>2 with different sound velocities, obtained by diagonalizing the quadratic 
form. In general however, the theory defined by 7i is not the sum of two independent 
conformal field theories. The massless excitations identified previously have to be thought 
of as a way to define the excitations of the whole theory, which is therefore not conformal 
invariant, since it does not have a well defined sound velocity. It would definitely be of 
interest to study the finite size spectrum of the lattice model to push this investigation 
further. 

The foregoing results generalize easily to arbitrary N, K. The key point is that the e 
excitations give rise to a massless U(l) (charge) degree of freedom, and that this excitation 
in turn feeds source terms on the I = 1 strings for both the (color) sl(N) and (flavor) sl(K) 
excitations. As a result, one always gets the SU(K) level one WZW model in the flavor 
sector, and a more complex theory in the color sector, that reduces to a mixture of the 
SU(N) level j WZW model and an SU(N) coset model in limit of small A. 



9. Impurities 

By a very general construction, it is possible to build an integrable impurity model by 
inserting a different representation, or a representation with a different spectral parameter, 
in the general quantum inverse scattering framework. This is the same trick that has often 
been used to study models with mixtures of representations, as well as models with spectral 
parameter heterogeneities |[4Q|| , ||26|| . In the context of impurities, the method was probably 
first used in |41]]. It was applied for instance to the t-J model with a four dimensional 
impurity (the si (2/1) case) in |42[ . 



Note that I am only discussing impurities in a periodic chain here; this is not the same 
(although results in the continuum limit are quite related) than having a chain with eg 
open boundaries, nor boundary impurities. For some recent results in that direction, see 



eg @ 



9.1. Fermionic impurities in a generic fermionic chain 

Let us first suppose that we insert in a chain based on fermionic representations with 
t > N — 1, a fermionic representation with a different value of the Dynkin parameter, 
and a shifted spectral parameter. The Bethe equations now look like (|2.2|) , except that 
for the /io root, the left hand side contains an additional term et>(no — T). The energy 



23 



takes the same form as before, and all the impurity term changes is the equations for the 
densities: equations ( |2.17|) and (|2.19| ) are unchanged, while (|2.16| ) contains an additional 
j^Gt'+i, p (no — T), L the size of the system. 

Since the equations for the densities appear in the thermodynamic Bethe equations 
only through their variations, the equations determining the ground state are unchanged: 
( |5.1|) , (|5.2|) and ( |5.3| ) still hold, together with the analysis of the previous sections. Pro- 
ceeding further to analyze the scattering of the excitations, equations (|6.5|),(p77|) and (|6.9| ) 
still hold, too. The only role of the impurity is that p^ and p^ in the left hand sides 
of ( |6.5| ) and ( |6.9[ ) have now a 1/L part, which follows from the solution of the equation 
generalizing ( |5.31| ): 

rQ 1 
p N (X) + K(X - p) PN (p)dp = G t+hN (X) + -G t , +1 , N (\ - T). (9.1) 

J — Q 

Note that the Fermi cut-off Q is not changed, since it follows from the condition e N (Q) = 0, 
and the latter is a TBA equation, independent of any impurity terms. The Fermi velocity 
has a 1/L correction that we neglect in the limit L large. The only effect of the impurity 



is thus to modify pN and pn in the equations (|6.5|) and ( |6.9|) by a term of order 1/L. 
Consider for instance equation (|0|). Since p~N still vanishes for rapidities smaller or equal 
to Q, the impurity term does not introduce any non trivial phase shift for the densities p p ; 
its only effect at leading order in the scaling limit Q — > oo is to renormalize the mass of 
the excitations (keeping their ratios constant) by a 1/L term which becomes negligible in 
the limit L large. The impurity does not introduce any flow in the renormalization group 
sense, and roughly corresponds to changing the length of the system by a finite amount 
(I'll refer to this, not quite correctly, as being "irrelevant") . The same conclusion holds 
for the u(l) sector. 

In contrast, suppose now that t' < N — 1. In that has been noticed before, 

the combination 

H pN G t '+i,N _ n G P NGt'+i,N 

<^t>+i, P 77 - l*t'+l, P 7—7; 

L + Unn L + LrNN 

does not vanish. Calling this combination If, pi it follows that for large L the right hand 
side of equation (|6.5| ) contains now two source terms, s p n * Pn\ and j^It' , p (p — T). 

In the case t' = 1, the impurity term is ^s P) at(/x — T): the equations exactly coincide 
with those of the exactly screened su(N) Kondo model with K channels ||44j| , ||45|| . although 
the bulk now appears massive. If however one concentrates on the massless limit of the 
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bulk degrees of freedom by letting the rapidity p — > oo, while at the same time also sending 
the impurity rapidity T — ► oo, the Kondo equations are then exactly recovered, 2nT/N 
being related in a simple way to the Kondo rapidity. 

In fact, it is easy to see that our Bethe equations coincide with those of the degenerate 
Anderson model when K = 1 - this has already been observed by Schlottmann in some 
cases [46]. Indeed, the equations for the degenerate Anderson model as written say in [44] 
are the same as the ones we are considering, with t' = l,T = |^,/xo = ^, except that the 
source term for the Anderson model is exp(ikjL) = exp(2iTfj,oL) instead of our [et(fio)] L - 
To match these too, it is enough to send t — > oo with t oc 1/r, and rescale the length of the 
system appropriately. Since the physical properties are independent of t for t > N — 1, we 
should indeed obtain the same results as for the degenerate Anderson model when t' = 1. 

Introducing the physical rapidity 9 = such that the dispersion relation of massless 
excitations is p = e oc e e , we find that 

m _l d smj^-^) _ld 

SN ~^ ~ 7^ sin( 4 + ^) = ~iTe {p) - (9 - 2) 

In the case of higher values of t', we find more complex reflexion matrices. For t' = 2 
for instance, one has (p — l)(p+ 1), etc. The meaning of these scattering matrices will be 
discussed further in J4"7| . 

To conclude, we see that, while the continuum limit without impurity was a Gross 
Neveu model with iV colors and K flavors, either the fermionic impurity is irrelevant if 
t' > N — 1, or it affects the su(N) sector in the same way as in a pure su(N) theory (as 
for the su(K) and u(l) sector, they are still unaffected). 

Non fermionic impurities in the fermionic chain could also be considered, but their 
effect is quite straightforward: they affect the SU (N) or SU (K) sectors of the Gross Neveu 
model as in pure SU(N) or SU(K) theories, giving rise to various N or K channel Kondo 
models (a solid state physics model corresponding to this situation was proposed in 



9.2. Fermionic impurities in a non fermionic chain 

Another example of interest is a fermionic impurity in a non fermionic chain. Consider 
for instance again the case N = K = 2 with a representation having Dynkin parameters 
a± = j in the bulk. The only difference with the analysis of section 8 is that a 1/L term 
appears in the right hand side of equation ( |8.3| ): the situation is thus very similar to the 
case of a fermionic representation in a fermionic chain: the impurity renormalizes the 
source terms for the two SU{2) systems by a term of order 1/L is rapidity independent: 
no flow is generated. This conclusion generalizes to any N, K, and seems to agree with 
the results of [g2j. 
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10. Conclusion 

In conclusion, it does not seem possible to observe any interesting "supersymmetric" 
properties in the continuum limit of integrable lattice models based on sl(N/K) superalge- 
bras . In fact, with the whole superalgebra symmetry, the continuum limit of these chains 
is not even relativistic, in sharp contrast with what happens in the case of ordinary alge- 
bras, where this continuum limit coincides with Wess Zumino models on the corresponding 
group. Interesting continuum limits can be obtained only when the superalgebra symme- 
try is broken. The case we have considered in details here leaves the sl(N) ® sl(K) (g> u(l) 
symmetry, and, in the continuum limit, gives various instances of color and flavor Gross 
Neveu models. The most "symmetric" case is obtained with fermionic representations with 
Dynkin parameter t > N — 1: in that case, one truly gets the iV colors, K flavors, Gross 
Neveu model, with the remarkable property that all the excitations have the same sound 
velocity indeed. Other cases lead to continuum field theories with less symmetry, typically 
involving mixtures of massive and massless excitations with different sound velocities. We 
have also considered impurity models, concluding that in that case too, nothing really new 
is observed, impurities either leading to irrelevant perturbations, or reproducing known 
Kondo models in the SU(N) or SU(K) sectors of the GN model. In particular, the con- 
tinuous parameter that is our disposal when using typical representations does not give 
rise to interesting tunable parameters in the field theory limit; in general, it simply affects 
the sound velocity, or the overall mass scale. 

It must be stressed that models based on the quantum deformations sl q (N/K) would 
be relativistic, even without the introduction of a chemical potential (this is especially 
easy to see in the case of si (1/1), where turning on the quantum group deformation is 
equivalent to adding up the chemical potential A). It is not known what the continuum 
limit of these models is in general, nor what happens to them as q — > 1; this is actually a 
very intriguing question, on which I hope to report soon. 

An issue that is somewhat related is what would happen for models based on, for 
instance, alternating fundamental representation and its conjugate. Unlike in the sl(N) 
case, the conjugate of the fundamental does not behave like another fundamental repre- 
sentation. While, for instance, the s/(3) model with an alternance of 3 and 3 is integrable, 
I do not know whether the s/(2/l) model is, nor what the Bethe equations would look like. 

1 Models based on osp(N/2M) seem more promising, since, according to [pt9| , the ones based 
on the fundamental representation are conformal invariant. 
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This problem deserves more study, as it seems related with important issues in disordered 
systems [|50] . 
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